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We investigate bouncing solutions in the framework of tiie non-singular gravity model of Brandenberger, 
Mukhanov and Sornborger. We show that a spatially flat universe filled with ordinary matter undergoing a 
phase of contraction reaches a stage of minimal expansion factor before bouncing in a regular way to reach 
the expanding phase. The expansion can be connected to the usual radiation- and matter-dominated epochs 
before reaching a final expanding de Sitter phase. In general relativity (GR), a bounce can only take place 
provided that the spatial sections are positively curved, a fact that has been shown to translate into a constraint 
on the characteristic duration of the bounce. In our model, on the other hand, a bounce can occur also in the 
absence of spatial curvature, which means that the timescale for the bounce can be made arbitrarily short or 
long. The implication is that constraints on the bounce characteristic time obtained in GR rely heavily on the 
assumed theory of gravity. Although the model we investigate is fourth order in the derivatives of the metric (and 
therefore unstable vis-a-vis the perturbations), this generic bounce dynamics should extend to string-motivated 
non singular models which can accommodate a spatially flat bounce. 

PACS numbers: 



I. INTRODUCTION 

Observations, especially those of the CMB by WMAP H], 
strongly suggest the occurrence of a primordial inflationary 
period [2 3 4|. Inflation not only provides an explanation for 
the homogeneity, flatness and horizon problems of the stan- 
dard hot big bang cosmology, but it also offers a consistent 
mechanism under which metric fluctuations are stretched be- 
yond the Hubble radius with a nearly invariant power spec- 
trum. However, in spite of its successes, the theory of infla- 
tion does not solve the problem of the initial singularity. Al- 
though the weak energy condition is likely to be violated in 
such models, Borde, Guth and Vilenkin have shown |5 1 that 
inflating spacetimes are in general geodesically incomplete to 
the past 16121. 

Mostly inspired by the string motivated pre-big bang sce- 
narios |8 9|, bouncing models ([TOlinillllllllEj, i.e., mod- 
els in which the universe undergoes a phase of contraction 
followed by expansion, have been proposed as alternatives to 
the inflationary paradigm. A bounce could solve the flatness 
problem of standard cosmology if the contracting phase lasted 
much longer than the expanding one, and it could also solve 
the homogeneity problem by making the past light cone large 
so that thermalization could occur. 

The problem of the initial singularity, however, is still 
generic. In the pre-big bang scenario the cosmological field 
equations exhibit a new symmetry, the scale factor duality, 
which maps the pre-big bang, contracting dilaton dominated 
era (in the Einstein frame) to the usual Friedmann-Robertson- 
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Walker cosmology (post-big bang phase) \15 \. Nonetheless, 
it has been shown that the two branches cannot be connected 
to each other smoothly |16|. This means that, while the pre- 
big bang era has a future singularity, the post-big bang phase 
emerges from a past singularity. 

The issue of how metric perturbations could be affected by 
the bounce has also been addressed, mostly in the framework 
of General Relativity (GR). The low energy approach repre- 
sented by GR is natural if we consider the fact that high energy 
corrections coming from, say, string theory, are negligible be- 
fore and after the bounce. Of course, this is no longer the 
case if we want to fully describe the bounce mechanism itself. 
This uncomfortable situation has led to the postulate that, in 
analogy to other short transitions in standard cosmology, such 
as pre-heating ||17|| or radiation to matter dominated epochs 
IfTSl . the time scale of the bounce is such as to permit fluc- 
tuations to evolve through it in a scale -invariant way llsl [T9l . 
In the GR framework, however, this assumption is far from 
being generic. On the contrary, it has been shown |20| that 
large wavelengths do suffer the influence of a such cosmo- 
logical transition. Hence the need to understand the way the 
characteristic bounce time is constrained (or not) by the field 
equations. 

In addition, because GR forbids the bounce to occur as long 
as the null energy condition holds, setting up such an evolu- 
tion for the scale factor can be a challenge. For instance, in a 
scenario with a single scalar field and positive spatial curva- 
ture lfT3l . the spectrum in the large scale limit exhibits A;-mode 
mixing. Models without spatial curvature but with a generic 
scalar field - or A;-essence field [21] - as matter content have 
also being built lEI . In that case, besides the fact that physi- 
cal observables are still affected by the bounce, those scenar- 
ios all lie in the phantom sector. Thus, in order to connect the 
fc-bounce to an expanding radiation era, a decay mechanism 
similar to pre-heating would be necessary. 
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In this work, we neither look for regularizations of the pre- 
big bang scenario nor for some other matter/curvature config- 
uration in the classical realm. Rather, we will adopt an alterna- 
tive approach, focused on a modified, higher order derivative 
gravity model proposed in 123, 24J. In this approach, an effec- 
tive action for gravity is constructed in such a way that all cur- 
vature invariants are limited. This is done via non-dynamical 
Lagrange's multipliers whose potentials ensure that the the- 
ory approaches Einstein's gravity at low curvature and that all 
solutions are well behaved at high curvature. 

In a subsequent paper |25 1, as an attempt to regulate singu- 
larities in the pre-big bang cosmology, Brandenberger, East- 
her and Maia studied a non singular dilaton cosmology in the 
framework of the model presented in 1231 1241. They found so- 
lutions corresponding to a contracting, dilaton-dominated uni- 
verse which evolves toward a bounce and emerges as a Fried- 
mann universe. Here we propose another solution: we will 
consider a homogeneous and isotropic universe filled with or- 
dinary radiation. The choice of curvature invariants will dic- 
tate the dynamics for the scale factor, which develops from 
a regular bounce to a Friedmann expanding universe, ending 
up with a quasi de Sitter expansion phase which could mimic 
the present acceleration of the universe. We also address the 
problem of the duration of the bounce in this model, and show 
that, as opposed to the classical treatment 1 13 1, it turns out to 
be completely unconstrained by the field equations. 

Extensions of GR involving higher powers of curvature in- 
variants are well justified, as calculations of one-loop diver- 
gences in quantum gravity generate terms proportional to R^, 
Rl^ and rI^^^ 126,^. As shown by Stelle in 1977 [28 1, al- 
though such actions can lead to a renormalizable theory, they 
all have a shortcoming, namely, the presence of ghosts - de- 
grees of freedom with negative kinetic energy [29 , 30 1. This 
fact makes the theory highly unstable, in the sense that the 
vacuum (empty) state can decay into a collection of both pos- 
itive and negative energy states. Also, at the classical level, 
one should expect that such instability will lead to growing 
gravitational perturbations carrying both positive and negative 
energy modes. This phenomenon is known in the literature as 
Ostrogradsky's instability |31, 32|. It is in close relationship 
with the fact that the Hamiltonian, due to the presence in the 
Lagrangian of derivative terms of order greater than one, is 
unbounded from below. 

In this work we will consider a simple quadratic theory, 
namely, a theory for which the higher order derivative terms 
appear only linearly. As a consequence, the background field 
equations will be of second order at the background level. 
This fact, however, does not suffice to prove stability, as the 
perturbations would still possess higher-order equations of 
motion. 

This paper is organized as follows. In the next section, we 
outline the model and write the full equations of motion, in 
terms of both the cosmic and conformal times. In section III 
we show how the universe can pass through a non-singular 
bounce and then connect to the usual expanding phases of 
standard cosmology. The final phase approaches asymptot- 
ically to an ever-expanding de Sitter period. At the end of 
section III we show how the time duration of the bounce can 



be made arbitrarily short. We conclude in section IV. 



n. THE MODEL 

Our goal is to investigate bouncing cosmological solutions 
in a non-singular, higher order derivative gravity proposed 
in Refs. [23, 24 [ In this model combinations of the Rie- 
mann tensor are introduced into the action via non-dynamical 
Lagrange multipliers. The potentials associated to each of 
the multipliers ensure that the theory approaches the Einstein 
limit at low spacetime curvature, and that the solution to the 
field equations is non-singular (typically, a de Sitter solution) 
at high curvature. More generally, for the gravitational action, 
we have 



1 r\ ^ 



gd'^x. 



(1) 



where G^^ - ^mn^^ is the Newton constant defining the Planck 
mass Wpi, ipi {i -{,■■■ ,N) are Lagrange multipliers depend- 
ing on space-time coordinates, and 7*'^ are functions of the 
curvature invariants 



(2) 



which we would like to limit. 

In order to understand how these limits are implemented, let 
us restrict ourselves to the case where a single field if\ = is 
present, and only one invariant, = 7, is limited. Variation 
of the action ([T]l with respect to (/? then provides the constraint 
equation 



di^ 



(3) 



At low curvature, we demand that be small and the theory 
to approach GR. Hence the action in this limit should be ap- 
proximated by 



(4) 



which implies, from Eq. ([3]l, that the potential should behave 

as 



lim V (if) ~ If' 



2 3 
+ (p + ■ 



(5) 



On the other hand, at high curvature, which we take to mean 
1^1 » 1, the correction term for the Einstein-Hilbert action 
becomes important. The potential must then be chosen in such 
a way that the solution of the field equations approaches the 
de Sitter solution, thereby effectively limiting the curvature. 
From Eq. ([3]l, we see that this requirement can be fulfilled 
provided that 



lim V(<p) = 2A, 



lim I (r,rI„...) = 0, 



(6) 
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FIG. 1: Scalar field potential V((^) (full line) and its derivative 
(dashed line), functions of the scalar field tp, as used for the modi- 
fied gravity field equations. 



where A is a constant. Thus, the action in the high curvature 
regime can be written as 



hm 5„ 



1 



2 A) V^d'^x. 



(7) 



We thus recover the Einstein-Hilbert action with a cosmolog- 
ical constant term coming from the potential. A detailed dis- 
cussion of this de Sitter limit can be found in |24|. 

Here we are mainly interested in the particular case where 
the universe, otherwise filled with a perfect fluid such as, e.g., 
dust or radiation, is originally undergoing a phase of contrac- 
tion where the spacetime curvature increases. Since the cur- 
vature is now limited, the Universe must reach a minimum 
radius, at which point it should bounce and subsequently ex- 
pand. 

Such a cosmology should be realizable with a conve- 
nient choice of curvature invariants. Assuming a Friedmann- 
Lemaitre-Robertson Walker (FLRW) metric of the form 



As' 



-dt^ + (f) 



dr^ 



1 -7Cr2 



+ r 



(d0^ + sin^ edip^) 



(8) 



we find that singularity-free second order equations can be 
obtained with the curvature invariant 



/ [r, R^yR"") = R- .y3(4^J^^!^'^^^^^^ = 12 (// -H H^) , (9) 

where the Hubble expansion rate isH = a /a, provided that the 
potential obeys the twin conditions that V(Q) = 0, V'(Q) = 
and V'i+oo) — > constant. 

One simple example of a potential that satisfies the condi- 
tions spelled above would be 



(10) 



and it is represented on Fig. [T] Inserting iMs potential into 
Eq. (|3]l, we obtain 



H + W 



1 dV 



12 d^ 2(1 -H 



,6^3/2 ■ 



(11) 



Eq. (Ill tells us that, at the bounce where // = 0, we can have 



H > provided we choose a positive cosmological constant 
A. Besides, we want H to reach a maximum value after the 
bounce, after which it should decrease again - we put this ex- 
tra ingredient in order to illustrate the connection between the 
bounce era and the usual cosmological setting of a decelerat- 
ing (H < 0) Universe. In the next section we will show that 
this feature is realized for a range of initial conditions. 
The complete gravitational action we choose is therefore 



1 



d\. 



(12) 



where 



-Cgrav - R + ifi 



R-^{ 



^R^.yRt"' 



R^) 



-V(ip). (13) 



Variation of the action ( [T2] i with respect to the metric then 
leads to the modified Einstein equations 

^"V - ^g^ivR- g^iv'^l + iphR^lV- (flR)-y-f,+ glIv(<Ph)fp 



+ 2<fIpR^pRpy + {ifilpR^yf + g^y {ipIpRP'), 



- 2 [pIpRf/^\^,^^ + ^-g,yV {if) = 87rG,r,„ (14) 

where we have defined Ir = dl/dR and Ip = dljd [R^.yR'"'); 
in Eq. ([T4]i, we also added a matter component whose stress- 
energy tensor is T^^y. These equations should be solved to- 
gether with Eq. ( [Tl] i. 

Eqs. (14 1 are, in general, of fourth order in the metric 
variables, and therefore subject to Ostrogradsky's instabil- 
ity. However, in the special case of the maximally symmetric 
FLRW metric ([8]l, and for the particular choice of the invariant 
/ of Eq. (|9|, we find that the second derivative of the scale fac- 



tor appears only linearly in the Lagrangian (13 i. As a result, 



the background equations of motion stemming from (14 1 are 
second order. Indeed, for the time-time component we obtain 
the modified Friedmann equation 



hUI + 4(p) + 2Hip + 



7C 



V SttG, 



(15) 



with p the matter energy density of the perfect fluid. For the 
spatial components the equations again lead to a generaliza- 
tion of the usual GR case, namely 



ip + 4H(p + {H+ -H^ (1 -I- 4(p) + 



•K 



V 



2fl2 12 



-4nG,p, (16) 



where p is the matter pressure. This last equation can also 



be derived from Eq. ( 15 1 and the matter stress-energy tensor 
conservation. 

A simple way of obtaining these equations is by setting 
goo - -N^i^), T being a general time coordinate. This trans- 
forms the gravitational action into 



SttG, 



d'*xL(N,<fi), 



(17) 
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where y]y — sin 6/ Vl - 'Kr^ is the square root of the de- 
terminant of the metric on the spatial sections, and 



•K 

7^ 



1 

N 



da 
fldr 



+ 2- 



dfl 
fldr 



d(p da 
d7^2^^ 



(18) 



1.0 



0.5 



0.0 



Varying Eq. ( 17 1 with respect to (fi and yields Eqs. (Ill and -Q 5 



( 15 I respectively once one fixes the lapse function to unity, 
i.e., once one identifies t with the cosmic time t. 

For later use it is convenient to rewrite the equations of mo- 
tion in terms of the conformal time rj, defined in the usual 
way as d?; = dt/a. Denoting by a prime the derivative with re- 
spect to conformal time, /' = df/drj, we rewrite the constraint 
equation as 

12 difi 

and the generalized Friedmann equation as 



(19) 



V 8nG, 



(20) 



These last two equations can also be derived from variations 

N 



of (17 1, setting t - t], i.e. N - a. Together with matter 
energy-momentum conservation 

p' + 3'H{p + p) ^0, (21) 



these equations can be combined to yield 

V AnG 

'H' {I + A^p) + if" + TH^p' ^ a^ 



■(p + 3p) 



(22) 



Eqs. ( 19l, (20 1 and (21 1 provide a self-contained description 
for the bounce, to which we shall now turn to. 



III. BOUNCING UNIVERSE 

Before studying the different cases of interest, let us use 
the cosmological constant A to define a more convenient time 
variable t = Va?. Rescaling the Hubble parameter, the scale 
factor and the energy density through 



H - Va/z(t), a - Aa and s(t) - 



(23) 



we obtain a set of dimensionless equations as 
dtp _ ^ 1 

d7 ~ m 



-h{l +4(p)- 



2a^h 2h ' 



and 



dh _ 1 d'V 
dr 12 dip 

where the potential is now 'V = V/A. 

The matter stress-energy conservation then reads 

de 
d7 



+ 3h (e + P)^ 0, 



(24) 



(25) 



(26) 



where the pressure p is modified in the same way as the en- 
ergy density through !P = 8nG^p /(3A). These rescalings also 
apply for the conformal time equations, with 'H unchanged. 




-1.0 



FIG. 2: Phase-space diagram in vacuum. The critical points are 
shown as filled circles, and are found at the intersections of the curves 
corresponding to dH/dr = and di/j/dr = (solid lines). The dashed 
lines show some particular solutions on this phase space. The arrows 
depict the arrow of time. 



A. Dynamics of the field equations 

We begin the section with a dynamical analysis of the vac- 
uum equations with zero spatial curvature (TC = 0), in the case 



of the potential of Eq. ( 10 1. We find that the system 



dh 



dr 2(1-H^6) 
dip 



,6^3/2 



1 



dr eh yflT^ 2 



--/i(l+4^), 



(27) 



(28) 



has six critical points (the number of crossings of the curves 
defined by d/i/dr = and di^/dr = 0) - see Fig. 1. One of 
these critical points, at ((p,h) ^ (-0.3,0.21), is an asymptot- 
ically stable point - an attractor. This means that trajectories 
come to an end (at f — > oo) at this point, which is hence a de 
Sitter attractor. 

The point {(p,h) - (-0.3,-0.21), on the other hand, is 
an unstable spiral point - where trajectories emerge from at 
f — » -oo. It corresponds to an anti-attractor, in this case, 
a contracting de Sitter solution. The other critical points, to 
the left and right of Fig. 1, namely, {(p,h) ^ (-1.25,-0.27), 
(-1.25,0.27), (1.35,-0.22) and (1.35,0.22), are all saddle 
points. Some particular solutions for this system are sketched 
in Fig. |2] 

The most important feature of this phase diagram, for us, 
is that the origin (tp - Q,H - 0) is not a critical point - it 
coiTesponds to the crossing of the two branches of the curves 
dH/dt - 0. In fact, trajectories coming from below (H < 0) 
will necessarily have to go through the origin in a finite in- 
terval of time. But the trajectories from below all cross the 
same point (the origin), which means that (at least for the au- 
tonomous system we considered thus far) this point is some- 
how singular. This means, in our case, that the trajectories 
passing through the origin will experiment a discontinuous 
change in the rate of change of ^ - in other words, one cannot 
predict the value of dip/dt for trajectories that emerge from the 
origin at H - Q. 
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To see that this is the case, we look for solutions near the 
origin ((p = Q, H = 0) of the phase plane. In this limit, the 
dynamical system can be approximated by 



dh 



= /!- — + 

d(p h 



0[h^,ip^,hip). 



(29) 



We can rewrite the above equation in terms of h^, and then 
find a first integral of the resulting differential equation. The 
final result is that 



+ 1+ 4(^(1 +2^), 



(30) 



where A is an arbitrary integration constant. This means that 
the 2 dimensional dynamical system is still Lipschitz con- 
tinuous at the origin (see, e.g., Ref. |33|, pg. 53), which 
indicates that the 2 dimensional trajectories can be seen as 
projections of trajectories in a higher-dimension dynamical 



phase space. Hence, Eq. ( 30 1 indicates that although the two 
branches (H > Q and H < Q) cannot be simply connected 
to each other in the case of the 2 dimensional phase space, 
any enlargement of that phase space will lift the problematic 
degeneracies. 

The saddle points divide the phase diagram in different re- 
gions through their separatrices (not shown in Fig. |2] although 
they can be inferred from the particular solutions.) In par- 
ticular, the separatrices in the third and fourth quadrants ex- 
plain why there are solutions which never reach the H - Q 
axis. This is a general feature, in the sense that it does not 
depend on the matter content: indeed, it follows directly from 
Eq. (Ill, which is independent of the matter content for the 
constrained theory we are considering. However, the evolu- 



tion of the right-hand side of Eq. (Ill being coupled to the 



matter content through Eq. (15 i, one expects a similar gen- 
eral behavior in the Hubble parameter for universes with any 
type of matter content, but with only minor differences in the 
details of the phase diagram. 

Hence, we conclude that a bounce in this constrained grav- 
ity theory does seem quite generic, but unfortunately the con- 
tracting and expanding branches cannot be connected in a 
continuous (unambiguous) way to each other, because the 
dynamical system is two-dimensional. However, as we will 
show next, the introduction of even a minuscule amount spa- 
tial curvature TC, or of any kind of matter, both opening up 
the phase space into a third dimension, is enough to regular- 
ize this bounce, so that every physical variable in the system 
remains finite and no singularity occurs. 



B. Spatial curvature 

By inspection of the phase diagram of Fig. |2] we can see 
that for every trajectory in the lower plane {h < 0) there is 
a similar trajectory in the upper plane. The only problem is 
that apparently all these particular solutions cross at the origin 
{ip - Q,h - 0), which, sadly for the 2 dimensional dynamical 
system we considered in the previous Section, cannot make 
sense. As it turns out, spatial curvature is one way to alleviate 
this problem. 



Indeed, with even a small amount of spatial curvature, tra- 
jectories coming up from the lower-right corner {(p > 0, 
dtp/dr < Oh < 0, d/i/dr > 0) emerge from the origin on 
the upper left comer {tp < 0, d^/dr < h > 0, dh/dr > 0), 
and solutions coming up from the lower-left corner emerge on 
the upper-right corner, so that the particular solutions to the 
equations of motion are indeed perfectly continuous. This can 
only happen because the presence of spatial curvature turns 
our original (flat) 2 dimensional phase space into a 3 dimen- 
sional phase space, where h and a are effectively independent 
dynamical variables. This means that the trajectories of Fig. 
|2] which cross = and h - can now do so freely, as the 
third dimension (the scale factor a, shown on Fig. [3]l opens 
the possibility that the paths do not cross - in other words, in 
the 3 dimensional phase space these trajectories actually do 
not intersect. 

To see this, consider the dynamical system composed of 



Eqs. (19 1 to (21 1. The problematic equation, which leads to 



the difficulties exposed in the previous section, is Eq. (20i. In 
fact, by taking the limit "TY — > at the putative bounce we can 
see that this equation reduces to 



1 (a^y 



(31) 



The rate of change of the scale factor is naturally zero at the 
bounce ("H — » 0), hence, if the derivative of the field above is 
finite, then ip ^ ips and a — > ob at the bounce. This means 
that V{(ps) - 6'7C/flg, which then implies, for the potential 



given in Eq. ( 10 1, that 



97C2 



1 



(32) 



Clearly, as TC — > O'', tps — > 0, which means that the physical 
system itself is perfectly well-behaved for any nonzero value 
ofTC. 

Hence, we conclude that it is possible to connect the two 
branches of the theory (the contracting and expanding phases), 
provided that there is any amount of curvature, however small. 
Since the limit TC — > is well-behaved, we can do without it 
entirely, and join the pieces of the trajectories between the two 
branches, constructing zero-curvature bounce models. 



It is also interesting to notice that the condition ( 32 1 implies 
that there is a maximum value for the spatial curvature 



Aai 



(33) 



which is in fact perfectly legitimate in our constrained theory, 
since we limit the spacetime curvature at all times - and that 
includes the time of the bounce itself. These considerations 
are illustrated on Fig. |3] showing the time evolution of the 
scale factor, the reduced Hubble rate and the scalar field, for 
three scenarios with different values of the spatial curvature. 



C. Constrained gravity with a matter fluid 

Having established that, qualitatively, spatial curvature 
does nothing to our model, let us continue with the analysis 
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FIG. 3: Scale factor air) (left panel), reduced Hubble parameter /i(t) (center) and scalar field (p{T) (right panel) as a function of time for almost 
vanishing curvature {"K = 10 full line) and small but non vanishing curvature (TC = 10"^ and "K = 10"^, respectively dashed and dotted 
lines). It is interesting to note that the vanishing curvature limit allows for a so-called emerging Universe model in which the scale factor 
is essentially constant for a long period of time after which expansion (in fact inflation) takes place spontaneously. Pushing the calculation 
backwards in time shows that the vanishing of the Hubble rate is only stable in the actual limit '7C — » 0: expansion begins after a finite amount 
of time, but the amount of past-time during which h = can only be made infinite provided = strictly; otherwise, the Minkowski state 
always originates from a contracting de Sitter phase. 



of the field equations in the flat case, but now in the presence 
of matter. We will show in what follows that since the in- 
troduction of matter enlarges the dimensionality of the phase 
space, it similarly allows the trajectories in the contracting and 
expanding branches to be joined - whatever the amount of 
matter present. 

If p is the matter energy density, then the system will ac- 
quire one more degree of freedom, that is, our dynamical 



equations are Eqs. (24i, (25 1 and (26 1, which take the form 



dh 
d7 

d^ 
d7 
dfi 
d7 



2(1 +ip'=) 



,6^3/2 ' 



+ 



1 CO-* 

= --Ml +4^) + . 

- -3h(l + w)e. 



(34) 
(35) 
(36) 



where, as usual, w = p/p - f/s is the equation of state for 
the fluid (w - for dust, w - ^ for radiation, and w = - j for 
a curvature-like matter.) 

Before proceeding with the detailed analysis of the phase 
space of this dynamical system, let us concentrate on a partic- 



ular class of solutions of the system ( 34 1-( 36 1 which gives us 



some information on how matter affects the bounce. For this 
purpose, let us return to the previous vacuum case, Eqs. 
and ( [28] l. We have seen that, in our model, not only /z — » at 
the bounce, but also that dh/dr and dip/dr vanish there. The 
evolution at the bounce is therefore dominated by the higher 
time derivatives - which can be tuned, through the initial con- 
ditions, to small values, making the duration of the bounce 
very large. This means that, depending on the initial condi- 
tions of the Universe, the spacetime can stay in this quasi- 
Minkowski (h = 0) state for a long period of cosmic time 
before undergoing a phase of acceleration with dh/dr > 0. 
However, we see that in the presence of matter the Universe 



starts in a de Sitter contracting phase and quickly evolves to- 
wards the bounce, emerging in a super-inflationary period be- 
fore reaching the ever-expanding de Sitter phase. This is ex- 
actly similar to what we found for the curvature case and ex- 
emplified in Fig. |4] 
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FIG. 4: Time evolution of the Hubble rate in vacuum (full line) and 
with matter (dashed). When a small quantity of radiation is added to 
the system [p^ad (0) = 10"^ in this example], the scale factor evolves 
from a contracting de Sitter phase to the bounce. The final state 
corresponds to an ever-expanding de Sitter period. The results are 
exactly similar when, instead of radiation, one considers dust or a 
curvature-like fluid. The three cases are shown in Fig. |5] It is 
worth noticing that the presence of the fluid dramatically reduces the 
bounce duration so that an emergent universe model in this case is 
very unlikely. 

Notice that, as mentioned before, with matter we do not 
have a phase plane anymore - the phase space has three di- 
mensions. The possible trajectories that can be achieved in 
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FIG. 5: Time evolution of the densities when the fluid consists of 
dust (w = 0, full line), radiation (w = ^, dashed line) and curvature 
(w = -|, dotted line) Hubble rate in vacuum (full line) and with 
matter (dashed). 




FIG. 6: Effect of including matter on the scalar field as a function of 
time. The full line represents the scalar field evolution in vacuum, 
while the dashed line includes a dust or radiation component (both 
curves are indistinguishable). 



we considered many different potentials, and, as a matter of 
example, we show on Fig. |7] the phase diagram obtained for 
the choice 



2A(p 



(37) 



which is represented on Fig. [8] 
1.0 
0.5 
S 0.0 
-0.5 



-1 







FIG. 7: Phase diagram for the potential V2 given by Eq. l |37[ ( dis- 
played on Fig. [S] 




FIG. 8; The potential V2(<f) (full line) and its derivative (dashed) used 
to derive the phase portrait of Fig. |7] 



D. Bounce Scale 



the projected phase space ((fi,h) are essentially the same as 
those obtained in the 2 dimensional case of Fig. (j2]l. 

Thus we have been able to set up the cosmology for a homo- 
geneous, spatially flat and isotropic universe filled uniquely 
with dust and radiation and that behaves as follows: the uni- 
verse starts contracting until it reaches a minimum radius, 
when it bounces and expands afterwards. The expansion pe- 
riod right after the bounce coiTesponds to a super-inflationary 
period (d/i/dr > 0). When the Hubble parameter reaches a 
certain maximum value (which depends on the details of the 
spacetime constraints), it starts to decrease, similarly as in the 
usual FLRW universe, until it reaches either a de Sitter attrac- 
tor or some other regime as its final state. 

The behaviors discussed above are generic and do not de- 



pend crucially on the choice (10 1 for the potential. To see this. 



In this subsection our aim is to calculate the characteristic 
time scale of the bounce described by the above model. More 
specifically, we will show that, in contrast to the general rela- 
tivistic description lfT3ll22l . the typical bounce time scale can 
be made arbitrarily small in this category of theories. 

Let us start with the following expansion ifTSi for the scale 
factor around the bounce 



a (ij) = flo 



1 + 



1 



+ 6 



(38) 



where we have set the bounce time as the origin of conformal 
time, i.e., ?7bounce = 0. In Eq. (38 1, 770 provides the bounce 
characteristic time scale. Similarly, any other quantity with 
as a subscript means that this quantity is to be evaluated when 
the bounce occurs, i.e. for 77 = 0. As explained in Ref. |13|, 
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the parameters 6 and ^ control the amplitude of the deviation 
from the de Sitter-like bounce, whose scale factor is given by 
a (77) = flo -y/l + tan^ {t]/t]o), thereby explaining the factor 5/24 
in front of the fourth term in Eq. ( 38 1. 

In a similar way as for the scale factor, we assume that the 
scalar field (firj), its potential V{ri) and the energy density 
p irj) all have a Taylor expansion about t] - 0, that is we set 



(fi (77) = tfO + ifil 

= Po+P\ 



\'7o/ \ml vw/ 



+ P2 



+ P3 — I + 



■ (39) 
(40) 
(41) 



where the parameters (fj, V, and p, are to be determined below. 

In order to determine the defining parameters of the scale 
factor, and hence the bounce scale itself, it is sufficient to ver- 
ify in what manner they are constrained by the field equations. 



In practice, we insert the expansion ( 38 1 into Eqs. ( 19 1, (20 1 



and (21 1 and identify terms order by order in 77/770 <s: 1. This 
leads to 



Pi 

P2 
P3 



= 0, 

3 

= --(1 + w)po, 
- -36(1 + w)p(). 



(42) 
(43) 
(44) 



which expresses the first terms in the energy density. Note 
at this point that a symmetric bounce, having 6-0, leads 
to an even behavior for the density, as expected. The first 
equality merely expresses energy conservation at the bounce 
conformal time. 

The cosmic time bounce duration fo can then be evaluated 

as 



2 ^ 2 2 _ 

'0 - '^0'?() ~ 



dV 



(45) 



where the last term is the derivative of the potential with re- 
spect to (fi evaluated at the bounce conformal time 77 = 0; this 
remains a free parameter, as can be seen from inspection of 
the following relations 



Vo 



V2 = 



— - IbnG^po, 



«() 
dV 
d^ 



(36(pi + (P2) , 



(46) 
(47) 



relating the parameters for the scalar field potential. Note at 
this point that for negative or vanishing curvature TC, the po- 
tential needs to be negative at the bounce, so that, as in GR, 
a simple massive scalar field cannot lead to a bouncing phase 
unless the Universe is closed. 

At the next orders, we obtain relationship between the 
higher derivatives of the scalar field at the bounce and the 



energy density at that point as well as the asymmetry of the 
bounce, namely 



n = -^(1 +4(^o) + 6^^ 



(P3 = -6(1 +4-(fo) - -(fii. 



7C 



- 47rG^p()(l + w) 



(48) 
(49) 



and finally one can identify the fourth order term in the scale 
factor expansion ^ as 



V3 + 24kG^6pq (l+w)-246- 
^ 26(1 +4-(fo)-ip3 



(50) 



An important conclusion that can be drawn from these sim- 
ple relations is that the bounce scale can be made as small as 
one wants. Clearly, we see that the larger the derivative of 
the potential at the bounce or, put in another way, the larger 
the constant VA which defines the dimension of V (ip) is, the 
smaller the bounce characteristic conformal time scale will be, 
for a given bounce length scale aq. This result can be consis- 



tently checked by numerical integration of system (20i: con 



sider the vacuum equations of the last subsection, Eqs. (27 1 
and ( [28] l. At the point (ip -Q,h - 0), we recall that the deriva- 
tive of the potential is exactly zero. But, according to the rela- 



tion ( 45 I, this implies that the physical time duration fo = ao'7o 
of the bounce goes to infinity or, put in another way, that there 
is no bounce at all. This confirms the results for vacuum ob- 
tained in the last section. 



IV. CONCLUSION 

Cosmological inflation now probably deserves to be in- 
cluded in the standard picture of the Universe. However, as 
inflation cannot be directly observed, it is desirable to find 
other, challenging models that could similarly solve the prob- 
lems of non inflationary cosmology while at the same time 
making different predictions for data sets - some of which 
may yet be observed as, e.g., non gaussianity or tensor modes 
of perturbations. Such a non-inflationary model could include 
a phase of contraction followed by a bounce, leading to the 
current expansion. 

Having a bouncing phase in the early Universe is a highly 
non-trivial demand in the framework of general relativity with 
well-behaved matter content. It usually requires a strictly pos- 
itive value for the spatial curvature, in contradiction with the 
current observational data. Therefore, unless a phase of in- 
flation is invoked after the bounce occurred, either the matter 
content or the gravity theory must be changed. 

The simplest way for a bounce to take place in a regular 
matter theory consists in demanding a positive spatial curva- 
ture to compensate for the positive energy density of matter, 
allowing the Hubble parameter to vanish at the bounce. But by 
doing so, it was found that the time duration of the bounce was 
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bounded from below so that an arbitrary short bounce could 
not take place. It was suggested that this was entirely due to 
the curvature, whose value was indeed crucial in limiting the 
bounce duration. 

In this paper our main purpose was to investigate a non sin- 
gular bounce in the framework of the modified gravity model 
of f24l. We found solutions for which the universe is de- 
scribed by a period of contraction followed by an expanding, 
super-inflationary phase before connecting to a more usual 
radiation- or matter-dominated epoch of the FLRW universe. 
The final period then emerges as an ever-expanding de Sitter 
universe. We have computed the duration of the bounce for 
that model, and we found that, whatever the bounce, irrespec- 
tively of whether or not it is symmetric, whether or not there 
is spatial curvature (independently of the sign of this curva- 
ture), and whether in vacuum or not, its typical duration is 
unconstrained. 

This arbitrariness in the duration of the bounce could 



have important consequences for the evolution of fluctuations 
through the bounce. However, although the background mod- 
els are ruled by second-order equations, perturbations will al- 
most surely show signs of the underlying instability of our 
higher-order gravity model, so at this stage it is not clear that 
perturbations in the framework of our models could be sensi- 
ble. 
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